Dissipative spin-boson model and Kondo effect in low dimensional quantum gases 
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We show that ultracold atoms confined in a quasi- Id trap 
can be used to study the dissipative spin-boson model with 
tunable tunneling splitting and strength of dissipation. It is 
shown that with realistic parameters this system allows to 
study the crossover from coherent Rabi dynamics to incoher- 
ent tunneling. 

An effective two state system (spin 1/2) coupled to a 
bath of harmonic oscillators is one of the most important 
models to describe the effect of dissipation in quantum 
mechanics [1,2]. The associated spin-Boson Hamiltonian 
was extensively studied for a wide range of problems, 
most notably the issue of decoherence for superpositions 
of macroscopically distinct states [1] or various condensed 
matter realizations like the tunneling of interstitials in 
solids which are coupled to phonons or to electron-hole 
excitations in a metal [5] . In spite of its apparent simplic- 
ity, the spin-Boson model exhibits a very rich behavior, 
which crucially depends on the low frequency weight of 
the effective environment spectrum. It ranges from a 
simple damped oscillation described by the phenomeno- 
logical Bloch equations well known e.g. from NMR to 
complete localization in the so called ohmic dissipation 
case [1]. The latter is in fact closely related to a gen- 
eralized Kondo-Hamiltonian, a paradigmatic model for 
quantum impurity problems. 

In most of these applications, the precise strength and 
frequency spectrum of the environment depends on mi- 
croscopic details and cannot be changed externally. This 
is the case e.g. in the context of macroscopic quantum co- 
herence, where the final observation of coherent superpo- 
sitions of counter-propagating currents in a SQUID [3,4] 
requires extremely low dissipation whose origin is still not 
completely understood on a microscopic level. Similarly, 
charged impurities in a metal are subject to ohmic dissi- 
pation through electron-hole excitations, however their 
strength is typically very small, depending on the precise 
scattering phase shifts at the Fermi energy. 

Taking advantage of recent experimental progress in 
cooling and trapping of low dimensional quantum gases 
[6,7], we show that cold atoms can be used to realize 
an ohmic spin-Boson model with a tunable strength of 
both the interaction and the effective tunneling ampli- 
tude. This allows to study the complete range of dynam- 
ics from coherent oscillations to incoherent dynamics in 
the Kondo regime and finally to localization. The oscil- 
lator bath in this case arises from the coupling to the 



low energy excitations of a Bose condensate, which pro- 
duces an ohmic spectrum in the one-dimensional case, 
where gapless quantum liquids exhibit universal low en- 
ergy properties. 

The spin-boson model describes a two-level system (a 
spin-1/2, characterized by the Pauli matrix operators a), 
interacting with a bath of harmonic oscillators (phonons) 
according to the following Hamiltonian: 
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where A is the "tunneling" amplitude (the Rabi fre- 
quency in the case of a "free" two-level system), hq and &J 
are the annihilation and creation operators of the phonon 
modes characterized by a momentum q and the disper- 
sion relation LOq = uq, where u is the velocity of sound. 
Formally eliminating the oscillator part of the Hamilto- 
nian, it can be shown that the reduced dynamics of the 
spin is completely determined by the effective density of 
states 
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for sufficiently small values of uj. Depending on the value 
of the exponent s there are several different regimes. The 
most interesting situations arises in the so called ohmic 
case s = 1. The latter is is equivalent to 1 /r^-Ising model, 
which is known to have a phase transition even in Id [10]. 

The spin-boson Hamiltonian (1) can be implemented 
in a setup involving a sample of trapped cold atoms. Con- 
sider the two distinct hyperfine states (let us call them 
the state a and b respectively) of the same atom confined 
within a quasi- Id trap (see the Fig. 1). The atoms of 
the two bosonic species are trapped by the two different 
external potentials, Va{x) and Vb{x), respectively. Both 
potentials are assumed to have a high frequency of the 
radial motion, so that the motion of the particles can be 
considered as quasi-ld (this implies that both the tem- 
perature T and the characteristic interparticle interaction 
do not exceed the characteristic frequency of the trans- 
verse confinement, ll!±). Most of the atoms are in the 
state a and form a dense quantum liquid (the number 
of atoms in the state a is large Na ^ 1). The poten- 
tial Vb is assumed to have the form of an optical lattice 
with very well separated and tightly confining potential 
wells. Both potentials should spatially overlap to allow 
the interaction between the atomic species. The opti- 
cal lattice should have fairly low filling: Ni, <C Na, so 
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that the wavcpackets of atoms b do not overlap with 
each other and hence have no cooperative interactions 
with the quantum Uquid. At the same time, Ni, should 
be sufficiently large, to facilitate (ideally non- destructive) 
measurement. 

The interatomic interaction is described by the three 
interaction constants gpp', with /?,/?' = a,b characteriz- 
ing the interaction between the atoms a, b and between 
the atoms a and b. All the interactions are assumed to be 
short-range. The values of coupling constants are then 
related to the 3d scattering amplitude in a fairly compli- 
cated way and may contain resonances. In the simplest 
case, when the scattering length, corresponding to the 
scattering of the components /3 and /3', appi, is much 

— 1/2 

smaller than the radial ground state size Z_l ~ ' 
(hereafter we use units such that fi = m = 1, where m is 
the mass of an atom), the effective Id "scattering length" 
is given by a simple relation: gpjsi = 2TTap0'/l'^. The 
value and even the sign of the interaction constants g/s/ji 
may also be changed by applying external magnetic field 
(Feshbach resonances, see [11] and references therein), or 
simply by changing the trap aspect ratio (i.e. by varying 
ui±, see [8]). 




FIG. 1. The suggested experimental setup. The 

Bose-liquid of atoms "a" (white circles) is confined in a shal- 
low quasi- Id trap (gray). The atoms "b" (black circles) are 
localized by a separate trap (dashed line). The transition 
between the atoms is performed by a laser light. 

The on-site repulsion between the atoms b is Ubb ~ 
gbb/h, where lb is the characteristic size of the ground 
state wavefunction ipb{x) localized within a single po- 
tential well of potential Vf,. The value of Ubb should 
be large compared with the characteristic energy of in- 
teraction within and with the Bose-liquid of atoms o: 
gaaP, \gabP\ < Ubb (but still Ubb ^ t^±)- Then, only one 
atom of the type b can be trapped on the same lattice site 
at the same time. Accordingly, any state of the atoms b 
can be modeled by the corresponding state of a spin- 1/2 
system, namely: the spin-up state corresponding to, say, 
the state with a single atom b, whereas the spin-down 



state corresponds to the situation with no atoms 6 on a 
certain site. 

To model the exchange interaction we suggest to put 

the atomic system into the far-detuned laser light, pro- 
viding a transition between the internal states a and b. 
The amplitude of the Rabi flip fl can, in principle, be 
spatially inhomogeneous. The part of the Hamiltonian 
describing the interaction of atoms b with the Bose-liquid 
can be presented in the form of the following spin Hamil- 
tonian: 

Hb = ^{l + a,) j dx7plpa + j dxn{^ai^b(T+ + h.c.), 

(2) 

where ^a(a;) is the annihilation operator valued- field of 
an atom a at the point x, Pa = ^a^a is the liquid den- 
sity operator. We note, that in order to ensure a sin- 
gle occupancy of the lattice site the interaction between 
the atoms b should be repulsive, gbb > and sufficiently 
large. At the same time, the interaction of the atoms a 
and b can have arbitrary sign (provided that the number 
of atoms b is sufficiently small). The first term of (2) 
contains a time-independent phonon displacement oper- 
ator (the identity part of 1 + az operator) , which can be 
always integrated out from the system's action and only 
shifts the effective Hamiltonian by an additive constant. 

Now we turn to the description of the atoms a. To 
implement the Hamiltonian (1) we need to exploit a sit- 
uation when the dispersion relation for the low energy 
excitations of the Bose liquid of atoms a is linear in mo- 
mentum, i.e. has no gap. In the long wavelength limit a 
Bose-gas is a Luttinger Liquid (LL), and is described by 
the Haldane hydrodynamics Hamiltonian [12]: 



j dxivjdcj? + VN{de - TTPaf), 



(3) 



where pa is the equilibrium liquid density, vj = irpa and 
Vn = k/ttp^ where k is the compressibility of the liquid. 
The fields (f> and 9 describe the superfluid velocity and 
the density fluctuations, respectively. 

The LL model (3) can be diagonalized by the following 
transformation: 
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9{x) = -iJ2 
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sign{q)e"'-{b, + bl^, (5) 



where L is the sample length. Accordingly, the Hamilto- 
nian (3) takes the form: 



Hll = Vs'^\q\blbg, 



(6) 
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where: Vs = (vnVjY/"^ = {n/pa)^^"^- As expected, the 
spectrum of the lowest energy excitations is characterized 
by a single sound velocity u = Vs (cf. Eq.(l)). 

The actual values of the phenomenological parameters 
Vs and K — (vj/vn)^^^ in the LL Hamiltonian (3) de- 
pend on Pa, the interaction gaa, and on the parameters 
of the external potential Va- The requirement of the ab- 
sence of a gap in the excitation spectrum sets a number 
of restrictions. First of all, the interparticle interaction 
has to be repulsive {gaa > 0). Second, the trapping po- 
tential should not allow the formation of the so-called 
Mott-insulator state. The Bose-field operator can be ex- 
pressed as 



de 
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exp(im0), (7) 



where the dimensionless proportionality coefficient de- 
pends on non-universal short interparticle separation 
properties of the Bose-liquid. On the contrary, the ex- 
pression for the density in terms of the field 11 does not 
have this ambiguity: 



Pa{x) = {pa +n(x)) ^ e-xp{i2m6{x)). 



(8) 
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In the case of a trapped Bose-gas, the Luttinger param- 
eter K depends on the interaction in the combination 
7 = 9aal Pa [12]. In particular, for weak interactions 
(7 — > 0) ^^(7) « 7r/7^/^ and thus can be very large. 
In the other limit, i.e. when the interaction is strong, 
7 ^ 00, iC(7) « 1 (Tonks gas limit, [9]). 

Eqs.(7) and (8) for the Bose-field operator and the 
particle density can be used to rewrite the interaction 
Hamiltonian (2). The role of the terms with various val- 
ues of TO is quite different. Since the phase B contains a 
quickly oscillating term irpaX, all the contribution con- 
taining exp(iTO0) with TO 7^ average out under the inte- 
gral sign in Eq.(2), provided that the localization length 
Zfc of the wavefunction ■0^ is sufRciently large: trpah ^ 1- 
Then, in the limit g — > we can put ga; « everywhere 
in Eqs.(4) and (5). Finally, keeping only the terms with 
TO = in Eqs.(7) and (8) we can derive the following 
representation of the Hamiltonian (2) : 
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+npl/'ll/\a+eM^m) + h.c), 



(9) 



where 
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and 17 17 (the unknown numerical coefficient is the 
combination of the two factors: explicit dependence of 



the result on the wavefunction ipb and the unknown factor 
in Eq.(7)). 

To prove the equivalence of the Hamiltonians (1) 
and (9) we follow the unitary transformation: H' = 
S-'^{Hll + Hi,)S with S = exp(cr2i0(O)/2) and identify 
M = Ws, A = n{palbY/'^, and 
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(2 A') 1/2 



(10) 



We note though, that the equivalence of the initial spin- 
boson Hamiltonian (1) and its quantum gases version (9), 
(6) can only be justified if for ql\j <^ 1. For larger values 
of q the interaction becomes more sophisticated, but also 
decreases as soon as qlb ^ 1. One can still use Eq.(9), 
assuming that the summation over q is restricted by the 



condition uj{q) 



< 



where the quantity u/lb plays 



the role of the high-frequency cutoff. 



a=l.l a=0 




FIG. 2. The atoms b population relaxation. In the coherent 
regime, at zero temperature, for a = 0.1 and a = 0. In the 
incoherent regime, at low temperature, for a — 1.1. 

The spin-boson Hamiltonian (1) can also describe an 
isolated spin-1/2 impurity dynamics in a 3d Fermi sys- 
tem (the so called Kondo model). A slightly generalized 
(anisotropic) version of the model can be formulated with 
the help of the following Hamiltonian: 



H 



E 



e{q)cl^Caq + J||5.s.(0) + J±^±sx(0), (11) 



where c^g is the annihilation of a fermion with (now 3d) 
momentum q and the spin projection cr, S is the spin 
of the impurity, s is the total spin of the electrons at a 
given point. The exchange couplings J^and J|| describe 
the interaction of the different spin components {S^s± = 
SxSx + SySy). The exact solution of the model is given 
in [13]. Mapping of the Hamiltonian (11) to the spin- 
boson model Hamiltonian (1) is described in every detail 
in [1,14], and is based on the fact that the spin density 
excitations of a free Fermi gas can serve as an oscillator 
bath. This leads to the identifications A ^ J±, u ~ vp 
with vp being the Fermi velocity of the Fermi gas, and 
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The behavior of the system (the character of the impu- 
rity spin dynamics) depends strongly on the sign of the 
coupUng J||: the situations with J|| < and J|| > 
are commonly referred to as the ferromagnetic and anti- 
ferromagnetic cases respectively. 

The analysis above shows, that the spin-boson model 
(1) can describe both the Kondo and our quasi-ld Bose 
gas with an impurity. On the other hand, this also means 
that the cold atoms version of the spin-boson model can, 
among other things, be used to model the Kondo Hamil- 
tonian and the related phenomena, provided that the val- 
ues of the independent parameters gab and K in (10) 
are selected accordingly, i.e. the couplings (10) and (12) 
match each other. Both systems are characterized by 
J{co) = 2au} in the u ^ limit, where 

In the Kondo language, a > 1 and a < 1 correspond to 
the antiferromagnetic and ferromagnetic cases, respec- 
tively. We note, that since both. K and Qat} are indepen- 
dent and practically unrestricted {gab can have arbitrary 
sign and practically any value, whereas I < K < oo), 
the parameter a can take arbitrary values, both small 
and large. For realistic values gab <^ Pa ct Ki 1/2K 
and is small. For example, for a potassium conden- 
sate of Na = lO"' particles the interaction parameter 
gab/na ~ 10~^, and K ~ 30, which means that a ~ 0.02 
and can be made larger by a Feshbach resonance. 

The dynamics of the "spin impurities" in our proposal 
can be observed by following the population dynamics 
of the atoms b in the presence of the laser light fl. In 
the absence of interaction between the "impurity" b and 
the LL of atoms a, the atom of impurity should undergo 
familiar undamped Rabi oscillations with the laser field. 
The actual damping is determined by the width of the 
involved atomic levels and will be neglected here. On the 
other hand, the interaction between the atoms a collec- 
tivize the Rabi oscillations of the "spin" and the collective 
degrees of the LL of atoms a. 

A particular interesting regime of the ohmic two-state 
system is the low T and a < 1/2 weak damping case. 
In this limit the occupation P{t) of, say, the state "6" 
exhibits damped Rabi oscillations [15] 

P{t) = cos (Art cos(r/)) exp(— Artsin(r/)), 

where rj = 7ra/(2(l - a)), = A(A/wc)"/(""^) (see 
Fig. 2). This result holds as long as aT ^ A^. At higher 
T the dynamics is incoherent (no oscillations should be 
visible). 

For 1/2 < a < 1 the system shows no oscillations and 
P{t) is a sum of exponentials [15]. This is the regime 
which is directly relevant for the Kondo model. The de- 
tailed description of the model dynamics is still unsolved. 



The situation becomes simpler again for a > 1, which 
is the ferromagnetic analogue of the Kondo problem. At 
T = the spin is localized, whereas at larger T there 
are incoherent transitions between the states at a rate 
~ T^"~^. This peculiar behavior was in fact observed in 
SQUIDs experiments [16] (see the Fig.2). 

In fact the suggested experimental scheme is not con- 
fined to a case of a single spin impurity in a quantum 
liquid. If the Nb is sufficiently large, the pseudo-spins 
can interact with "host" liquid in a collective manner. 
This situation is described by Anderson model and is 
characteristic to a number of important solid state sys- 
tems, like heavy fermionic compounds (see e.g. [13]). The 
spin-boson model is also a paradigmatic model for qubits 
damped by an environment. Understanding how such a 
system reacts under f,unahle conditions, is obviously im- 
portant for all those interested in quantum computation. 
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